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ABSTRACT 

A new, transient ultraluminous X-ray source (ULX) was recently discovered by 
Chandra in M31 with a luminosity at ~ 5x10'^^ erg s~^. Here we analyse a series 
of five subsequent XMM-Newton observations. These show a steady decline in X-ray 
luminosity over 1.5 months, from 1.8 — 0.6 x 10^^ erg s~^, giving an observed e-fold 
timescale of ^ 40 days. This is similar to the decay timescales seen in multiple soft 
X-ray transients in our own Galaxy, supporting the interpretation of this ULX as 
a stellar mass black hole in a low-mass X-ray binary (LMXB), accreting at super 
Eddington rates. This is further supported by the lack of detection of an 0/B star 
in quiescence and the spectral behaviour of the XMM-Newton data being dominated 
by a disc-like component rather than the power-law expected from a sub-Eddington 
intermediate-mass black hole. 

These data give the best sequence of high Eddington fraction spectra ever assem- 
bled due to the combination of low absorption column to M31 and well calibrated 
bandpass down to 0.3 keV of XMM-Newton in full frame mode. The spectra can be 
roughly described by our best current disc model, bhspec, assuming a 10 Mq black 
hole with best fit spin ~ 0.4, declining from L/LEdd = 0.75-0.25. However, the data 
are better described by a two component model, where the disc emission is signifi- 
cantly affected by advection, and with an additional low temperature Comptonisation 
component at high energies which becomes more important at high luminosities. This 
could simply indicate the limitations of our current disc models, though changes in the 
energy-dependent variability also weakly supports a two component interpretation of 
the data. 

Irrespective of the detailed interpretation of the spectral properties, these data 
support the presence of accretion onto a stellar mass black hole in a LMXB accreting 
in the Eddington-regime. This allows an unambiguous connection of this object, and, 
by extension, similar low luminosity ULXs, to 'standard' X-ray binaries. 

Key words: accretion, accretion discs - X-rays: binaries, black hole 



1 INTRODUCTION 

The broad population of X-ray sources in our own Galaxy 
includes, amongst others, X-ray binaries (XRBs) believed to 
be powered by accretion of material onto a compact object, 
usually a neutron star or black hole (BH). These dominate 
the X-ray sky of the Galaxy and are important objects to 
study as they provide high quality data, allowing an insight 
into the nature of accretion. In the case of low-mass com- 
panion systems (LMXBs) the X-ray emission characteristics 
are well described by the presence of an optically thick, ge- 
ometrically thin accretion disc (Shakura & Sunyaev 1973) 



and/or Comptonisation of the disc emission in a hot plasma 
(Titarchuk 1994). However, as the luminosity of these sys- 
tems approaches their Eddington limit our understanding 
breaks down. Proposed models have included the produc- 
tion of large-scale winds driven by the large radiation pres- 
sure (Poutanen et al. 2007; Kajava & Poutanen 2009; King 
2004) whilst others favour the disc structure becoming ra- 
diatively inefficient and instead dominated by optically thick 
advection processes (Abramowicz et al. 1988; Mineshige et 
al. 2000). Although a small number of XRBs have been 
inferred to accrete at these high rates (V404 Cyg: Zycki, 
Done & Smith 1999; GRS 19154-105: McClintock & Remil- 



© 2011 RAS 



2 M.Middleton et al. 




Figure 1. Left: XMM-Newton EPIC PN image of M31 ULX-1 with the 35" extraction region overlaid. Right: Chandra HRC-I image 
of the source taken eleven days earlier. Through performing relative astrometry between the Chandra and XMM-Newton fields, we are 
confident that these are images of the same source. 



lard 2006; V4641: Revnivtsev et al. 2002; Cir X-1: Done & 
Gierliriski 2003), the large intervening column of material 
within our own Galaxy prevents a detailed analysis of the 
behaviour of the disc. We cannot even look to the active 
centre of distant galaxies (AGN) for answers, as mass scal- 
ing places the disc emission firmly in the unobservable UV 
(although see Czerny et al. 2011; Done et al. 2011). 

The proximity of our nearest neighbour galaxy, M31, 
has allowed the latest generation of observatories to detect 
(and often resolve) the X-ray bright menagerie of objects it 
plays host to. The range of objects in M31 is analogous to 
the population of X-ray sources we see in our own Galaxy 
(e.g. Pietsch et al. 2009). However, unlike our Galaxy, obser- 
vations of the M31 population is not hindered by the large 
absorbing column of the Galactic disc (see Warwick et al. 
2011) and the distance to the source has a much smaller 
relative uncertainty, improving the constraints on the X-ray 
luminosity. 

One remarkable object in M31 has been seen to emit 
at very high luminosities based on the estimated distance 
of 0.7-0.8 Mpc (see Vilardell et al. 2010; Tanaka et al. 2010 
etc.). CXOM31 J004253.1+411422 (M31 ULX-1 hereafter) 
was first detected by the high resolution camera (HRC-I) on- 
board NASA's Chandra mission at a count rate of 3.28 ± 
0.04 ct (Henze et al. 2009, ATEL #2356). The follow-up 
Stui/t observation (OBSID 29762) placed a rough limit on the 
flux of the object in the Chandra observation at ~5x 10"^^erg 
s~^. Due to its luminosity (> 10^® erg s~^) and distance from 
the centre of M31 (> 2'), this is a bona fide ultra-luminous 
X-ray source (ULX, see Roberts 2007). A sequence of several 
XMM-Newton observations, spaced at 7-10 day intervals, fol- 
lowed this bright detection and have captured the emission 
characteristics of the source over a period of ~1 month (see 
also Kaur et al. 2011). 

Assuming that the spectral behaviour is similar to that 
of other low-luminosity ULXs (e.g. Middleton, Sutton & 
Roberts 2011; Gladstone, Roberts & Done 2009; Weng et al. 
2009; Dubus, Charles & Long 2004; La Parola et al. 2003; 



Parmar et al. 2001; Takano et al. 1994) then we expect there 
to be a significant disc blackbody component to the emis- 
sion. Whilst ULXs are generally persistent in nature (al- 
though see e.g. Ghosh et al. 2006, Sivakoff et al. 2008), this 
object is clearly transient and, on this basis, has much more 
in common with the larger population of Galactic XRBs. 
These observations may therefore provide an unobscured de- 
scription of accretion through a disc at high mass accretion 
rates. In this paper we will describe the temporal and spec- 
tral properties obtained from the high quality XMM-Newton 
data as the source dims from its peak observed brightness. 



2 DATA ANALYSIS 

In total, M31 ULX-1 has been observed 7 times in outburst; 
initially it was detected using the Chandra HRC-I, followed 
by a Swift target of opportunity (ToO) and then 5 XMM- 
Newton observations. The details of each of these are given 
in Table 1. To confirm that we are viewing the same object in 
all the observations we extracted the HRC-I data of the field 
when the source is bright (Fig. 1) and obtained the positions 
of all the point sources using CELLDETECT in CIAO v 4.1. 
Wc then used edetect in SAS vlO to obtain the analogous 
source positions in the EPIC PN field from the observation 
chronologically closest to the Chandra observation (OBSID 
0600660201: Fig. 1). The positions of several bright sources 
common to both fields were then used to constrain their 
relative astrometry. The position of the ULX (as determined 
by the high resolution HRC-I of Chandra) at R.A.(J2000) 
= 00:42:53.15, Dec.(J2000) = +41:14:22.9 (with a typical 
error of ~0.5") was consistent (within errors) in both fields, 
confirming it as a single transient object viewed at different 
epochs. 

For each of the XMM-Newton detections we extracted 
the event files using SAS vlO and filtered the products for 
standard patterns (< = 12 for MOS, <= 4 for PN) and flags 
(=0 for spectral and timing products). From these we ex- 
tracted spectral and timing properties from circular, 35" 
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Table 1. Observations of M31 ULX-1 



Instrument 


OBSID 


obs. date 


useful exposure (ks) 


/x 


count rate 








(xlO~^^ erg cm~2 s~^) 


(ct s-l) 


Chandra HRC-I 


10886 


17-12-2009 


19.0 


6.40* 


3.3 


Swift XRT 


29762 


22-12-2009 


4.4 


2.40* 


0.7 


XMM-Newton EPIC 


0600660201 


28-12-2009 


14.5 


2.52 


4.2 


XMM-Newton EPIC 


0600660301 


07-01-2010 


13.3 


1.87 


4.6 


XMM-Newton EPIC 


0600660401 


15-01-2010 


6.4 


1.37 


3.6 


XMM-Newton EPIC 


0600660501 


25-01-2010 


10.0 


1.14 


2.9 


XMM-Newton EPIC 


0600660601 


02-02-2010 


9.7 


0.86 


2.4 



Notes: Observations of M31 ULX-1 by various missions over a period of ~ 1.5 months. The useful exposure is the 
length of the observation following removal of background flares in the Swift XRT and XMM-Newton EPIC cameras 
(PN value given) respectively, /x is the unabsorbed, integrated flux between 0.2-10 keV (from applying a simple 
continuum model to the data) and count rate is given for the PN in the XMM-Newton observations. Fluxes denoted 
by an asterisk (*) are taJjen from the literature values (Henze et al. 2009). 



radius source and background regions from the same chip 
avoiding other sources in the field. We also filtered for hard 

(10-15 keV) background flaring intervals based on data from 
the entire chip, leaving the useful exposures given in Table 
1. 

In the first observation there was partial pile-up in the 
PN which we corrected for by removal of a 5" radius region 
around the ccntroid of the PSF, chocking that the observed 
patterns were within the error bounds of the model patterns 
using the tool epatplot. We found significant pile-up in the 
the first two M0S2 observations and, due to requiring large 
centroids for removal, we instead ignore these datasets. 

Wc extracted spectral products using xselect, obtain- 
ing background-subtracted spectra and responses. We fit 
these using comptt to phenomenologically describe the con- 
tinuum, with a neutral absorption column^. This describes 
the spectra well (x^ = 2809/2789) for electron temperatures 
ranging from 0.94 0.78 kcV, i.e. all spectra have a rollover 
well below 10 keV. We use this model to derive the unab- 
sorbed 0.2-10 keV fluxes reported in Table 1. 

We plot the unabsorbed flux as a function of time in 
Fig 2. There is a clear exponential decay, similar to that 
often seen in X-ray novae which can be modelled as fit) = 
/o xe"*'''^ where / is the flux after t days, /o is the initial flux 
and r is the e-fold time. Given that the Chandra flux is only 
a rough estimate (Henze et al. 2009; ATEL #2356) we fit 
the full dataset (blue dashed line in Fig. 4) and the dataset 
without this point (red dot-dashed line) with a linear model 
in natural log space. These give e-fold times between 29 and 
39 days respectively, consistent with that observed for X-ray 
novae in our own Galaxy (Shahbaz et al. 1998). 



3 X-RAY SPECTRA 

The XMM-Newton observations range from Lx=1.8- 
0.6x10^® ergs s~^, so only the first few datasets meet the 
ULX criteria of >10'^^ ergs s~^. These luminosities and 



^ with a lower limit based on the Galactic line-of-sight to M31 of 
~6.7xl020cm-2 (Dickey & Lockman 1990) 



heavily curved spectra (see previous section) clearly re- 
late to standard stellar mass BHs accreting at close to 
the Eddington limit. The alternative explanation of an 
intermediate-mass BH (IMBH: Colbert & Mushotzky 1999), 
with 10^~^ Mq, would instead have L/LEdd <S 0.1 where 
the spectrum would be expected to be dominated by a hard 
power-law extending to ~ 100 keV. This is plainly inconsis- 
tent with the observed rollover in the data below ~10 keV. 
However, the Chandra discovery luminosity is ~ 5x lO^''^ ergs 
s~^, so this source is clearly connected to the ULX regime. 

Ruling out the presence of an IMBH means that the 
compact object driving the accretion is most likely a stel- 
lar mass BH. Assuming that the BH population is similar to 
that of our own Galaxy then we could be observing accretion 
onto a BH of mass ~5-20 M© at rates of L/LEdd ~ 0.2 — 2. 
XRBs accreting at similar rates have spectra which are dom- 
inated by a hot accretion disc, showing temperatures of 
~ 1 kcV for luminosities up to Eddington (e.g. McClintock 
& Remillard 2006). In these disc dominated states, the maxi- 
mum disc temperature and total luminosity change together 
such that L oc T*, indicating a constant inner omitting area, 
as predicted by the existence of the innermost stable circu- 
lar orbit in general relativity (see e.g. Done, Gierlinski & 
Kubota 2007). 

3.1 Disc models 

The most accurate available model for disc emission is bh- 
SPEC (Davis ct al. 2005) which self-consistently determines 
the spectrum at each radius by solving the vertical radiative 
transfer, and propagates this to the observer through the full 
general rclativistic spacetime. These spectra are integrated 
over the entire disc, assuming the stress free inner bound- 
ary condition, with innermost radius fixed to the innermost 
stable circular orbit. 

We fit this together with neutral absorption^ to all the 
PN data across all of the observations simultaneously (wo do 
not use the MOS datasets due to pile-up issues and the low 
effective area at high energies relative to the PN), fixing the 
inner radius of the disc (set by the BH spin) and inclination 
to be the same in all the datasets. As this source looks similar 
to many soft X-ray transients in our own Galaxy, we first 
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Table 2. Best fitting spectral parameters. 



Observation 


1 


2 3 


4 


5 


TBABS*BHSPEC 


log L/L^^d 


n 1,q +0.016 

-U.lOO _o.o26 


riH = 0.067xl022cm-2 
cos(i°) = 0.87 

a. = 0.36 t°o.{'i 
BH mass (Mg) = 10* 
-0.283 IHl? -0.418 t° °\t 


n t;91 +0.005 


-0 612 +° °°^ 
U.DIZ _o.oo6 



(d.o.f.) 3169.4 (2801) 
Null P 1x10-6 



log L/Le 



n ')'i7 +0.015 



nn < 0.068xl022cm-2 
cos(i°) = 0.86* 
a* >0.76 
BH mass (Mq) = 14.9 tol 



-0.403 



^0.016 
-0.005 



-0.537 



+0.016 
■0.006 



-0.640 



+ 0.016 
-0.005 



n -,oi +0.016 



(d.o.f.) .3078. .5 (2801) 



Null r 2x10 



log L/ LEdd 



-0.381 



+0.007 
-0.003 



riH 



= 0.067xl022cm- 
cos(i°) = 0.5* 



: 0.76 



+0.01 



BH mass (M©) >29.6 

hO.006 n cc/i +0.008 



-0.530 



-0.664 



n TCC +0.007 

-0.766 Io,oo5 



n otry +0.008 

-u.oor -0.005 



(d.o.f.) 3009.5(2801) 
Null 1' o.oo:-! 



TBABS* (DISKPBB+COMPTT) 



nn 0.104 +°-°°\ xl022cm- 



(keV) 



P 

norm 

^'^compjSccd (kcV) 



comp 



(keV) 



+0.08 
-0.08 
+0.03 
-0.01 

7.59 

fO.lO 
-0.13 
+0.30 
-0.10 

11.33 ttil 



0.44 
0.60 



0.55 
0.98 



+0.01 
-0.01 
+0.03 
-0.02 

3.69 

f0.04 
-0.04 
+0.03 
-0.02 

16.39 +tf. 



0.51 
0.60 



0.96 
0.49 



0.67 
0.60 



+0.03 
-0.03 
+0.02 
-0.01 

1.14 

61 

'^■"^ -0.16 
n ai +0.06 

>19.00 



0.80 
0.56 



+0.03 
-0.03 
+0.01 
-0.01 

0.37 

<1.26 
<11.73 
(25.50) 



0.92 
0.55 



+ 0.03 
-0.03 
+0.01 
-0.01 

0.18 



(d.o.f.) 2752.2 (2777) 
Null P 0.73 



Notes: bhspec model: best-fitting values of the Eddington fraction are given for spectral fitting, with 
fixed parameters denoted by *. We firstly fit the data with the mass fixed in the model at 10 Mq (to 
be consistent with the BH population of our own Galaxy) with the inclination left free to vary, and 
subsequently with the inclination fixed at 30 and 60 degrees with BH mass free to vary. DISKPBD model; 
best fitting parameters arc given for the 5 observations with the 5*'^ having no constrained Compton 
component in the spectrum. Where upper or lower limits are given, the lower (0.01 kcV for the seed 
photon or electron temperature) or upper (200 for the optical depth, r, 0.8 for the spin, a, and 30 Mq 
for the BH mass) hard boundaries for the model parameter have been reached. In the case of the optical 
depth for observation 4, the value is given in brackets as it is unconstrained. In all cases where there 
are two components, the degenerate nature of the spectral fitting requires one best-fitting component to 
be held in position whilst 90% error bounds are placed on the other. This is an unavoidable result of 
modelling such broad continua. 
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Figure 2. Unabsorbed 0.2-10 keV fluxes of M31 ULX-1 from observations taken with Chandra, Swift and XMM-Newton (see Table 
1 for details of each of these). We fit the data with a In-linearly decreasing trend and find that, by including the Chandra detection 
(blue dashed line), we obtain an e-fold time of 29 days and excluding it (red dot-dashed line), an e-fold time of 39 days. These derive 
from spline fitting in QDP and have variances of 0.36 and 0.03 with and without the Chandra flux point respectively. Both of these are 
consistent with the observed e-fold times of X-ray novae (e.g. Mineshige, Yamasaki & Ishizaka 1993). 



assume a 10 M© BH. This gives a poor fit to the data (x^ = 
3169/2801) while returning a reasonable value for the spin 
and inclination (see Table 2). However, we obtain a better 
fit by allowing the BH mass to be free. We tabulate results 
for i = 30° and 60°, with the latter giving the best fit (x^ 
= 3009/2801) and a BH mass pegged at the upper limit of 
30 Mq. We show the data and residuals deconvolved with 
this best fit model in Figure 3. 



The residuals indicate that the spectrum is somewhat 
broader than the disc models (resulting in the poor fit qual- 
ity). This is commonly seen in XMM-Newton spectra of 
disc dominated XRBs (Kubota et al. 2010; Kolehmainen 
et al. 2011), but the level of residuals here is about twice 
as large. This could indicate that the wider bandpass used 
here (down to 0.3 keV rather than the 0.7 keV lower limit 
from the other CCD observations of disc dominated spectra) 
enhances the differences already seen between our best disc 
models and the data. Alternatively, the residuals could indi- 
cate that there is an additional component in the spectrum. 
XRBs commonly show a power-law tail to high energies, so 
we allow for this using the SIMPL model for Comptonisation 
(Steiner et al. 2009), fixing the power-law photon index to 
2.2 as is commonly seen in the high/soft state (Kubota, Mak- 
ishima & Ebisawa 2001). We include this additional compo- 
nent in our best fit model, but the improvement in is not 
significant. The scattered fraction is always < 5%, and the 
best fit disc parameters are not significantly different. 



3.2 Disc plus low temperature Comptonisation 

As the disc approaches the Eddington limit, optically thick 
advection should become important, where photons are 
swept along radially with the flow rather than escaping ver- 
tically (slim discs: Abramowicz et al. 1988). These photons 
can then be released in the plunging region where the flow 
becomes optically thin as it accelerates towards the event 
horizon (S^dowski 2009). Another process which can be- 
come important is a radiatively driven wind (Poutanen et al. 
2007; Ohsuga & Mineshige 2011). At very high luminosities 
we expect the wind to launch at large radii in the disc, but 
at lower luminosities we may expect the wind to be launched 
only from the innermost radii of the disc where the radia- 
tion pressure is strongest. These processes are not mutually 
exclusive, and it seems most likely that both advection and 
wind-driving affect the spectra at Eddington and beyond 
(Ohsuga 2007; Ohsuga et al. 2009). 

We have previously suggested that a wind-dominated 
model can explain both the spectral and timing character- 
istics of other ULXs at similar luminosities (Gladstone et 
al. 2009; Middleton et al. 2011b). The evolution of such sys- 
tems is difficult to directly observe as most are persistent. 
However, in this case we see a clear drop in luminosity al- 
lowing the predictions of the model to be tested. In par- 
ticular we expect that, as the mass accretion rate drops, 
the outer photospheric radius of the wind moves inwards 
as the radiation pressure decreases. The cool disc beyond 
the wind then has a smaller inner radius and should appear 
hotter. We describe this in xspec by a model comprising 
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a slim disc with the index of the temperature profile, p, a 
free parameter, together with thermal Comptonisation (with 
the seed temperature allowed to be free rather than being 
fixed to that of the disc) and neutral absorption^ (in xspeC: 
tbabs*(diskpbb+comptt)). Although there is no real ad- 
vantage to fitting the data simultaneously, we do so to allow 
a direct comparison with the disc-only model (bhspec). 

The best-fitting model and associated residuals are 
shown together in Figure 3 (right hand panel) with the 
model components highlighted in Figure 4 and the par 
rameters given in Table 2. We obtain an overall signifi- 
cant improvement in fit quality compared to the single disc 
model (Ax^ of 257 for 24 extra d.o.f.)- The best-fitting 
model parameters show that the temperature of the disc 
increases with decreasing luminosity, while its radius de- 
creases. The corresponding Compton component peak tem- 
perature (which is a function of both seed photon tempera- 
ture and electron temperature) also increases, but the frac- 
tion of the total luminosity carried by this drops with lu- 
minosity. This is broadly consistent with the predicted be- 
haviour of a wind/pliotosplicrc launched from smaller radii 
as the radiation pressure drops (Middleton ct al. 2011b). A 
possible issue arises where the disc at the lowest observed 
luminosity appears heavily advection dominated which is in- 
consistent with the inferred sub-Eddington mass accretion 
rate. Part of this is due to the lack of rolativistic smearing in 
this model which makes the disc emission artificially narrow 
(Kolehmainen et al. 2011; Kubota et al. 2010). 

We can obtain a rough estimate for the inner radius of 
the disc from the normalisation of the DISKPBB model. This 
drops from an apparent radius of 300 to 47 km assuming an 
inclination of 60°, which corresponds to a true inner radius 
of ~400 to 62 km after correcting for colour temperature and 
stress free inner boundary condition (Kubota et al. 2001). 
The lowest number is similar to that derived from the BH- 
SPEC fits for a 10 M0 BH i.e. 71 km for a = 0.36, whereas 
the larger radii seen for the disc at higher luminosities in this 
model imply that the wind photosphere dominates down to 
26Rg. This is somewhat larger than expected given that a 
10 Mq BH is still sub-Eddington at this point, so should 
not yet power a strong wind. This would argue for a smaller 
mass BH (e.g. XTE J1650-500, Orosz et al. 2004), in order 
that the first few observations are super-Eddington. 



4 SHORT TIMESCALE VARIABILITY 

The spectra alone are somewhat degenerate. As wc have 
shown, a description using only a single disc component 
(with fixed inner radius as the luminosity declines) is sta- 
tistically poor. However, the best available models are not 
sufficiently well calibrated with the expectation of residuals 
when compared to the data. An alternative (and statisti- 
cally better) description is a two component model, with an 
advective slim disc at low energies, together with an inner- 
disc wind photosphere at high energies which decreases in 
radius as the luminosity declines. We now use the additional 
information contained in the lightcurves to try to discrimi- 
nate between these two very different interpretations of the 
spectra. 

We can examine the variability as a function of energy 
within a given observation in several ways. The simplest of 



these is to obtain the variance of a light curve in a given 
energy band, subtract the expected variance of the Pois- 
son (white) noise and normalise by the moan count rate 
to get the fractional variability (excess variance, or rms, 
see Edelson et al. 2002, Vaughan et al. 2003). We split the 
bandpass into 4 separate energy bins, from 0.3-0.7, 0.7-1.4, 
1.4-3.0 and 3.0-10.0 keV and create lightcurves in each en- 
ergy range, bimiod on 300 s (to provide adequate statistics 
for such an analysis) from coadding the PN and MOS data 
(we note that, due to the lower quality statistics, the small 
amount of pile-up in the first two M0S2 datasets adds only 
a small constant offset to variability at hard energies and 
so will have negligible effect upon the result). These showed 
some weak evidence for variability. We then carefully iter- 
ated the binning across the energy bandpass whilst ensuring 
adequate statistics remained in each bin under investigation 
and were able to find at least one energy range for each ob- 
servation in which variability was detected at more than 3a. 
However, the upper limits on variability in the remaining 
energy bins are large. 

We can improve the statistics by using the energy band 
in which the variability is detected £is a reference lightcurve 
to calculate the fractional covariance i.e. the amount of vari- 
ability in the other energy bands which is correlated with the 
detected variability (see Wilkinson & Uttley 2009). The re- 
sulting energy-dependent fractional covariance spectra are 
shown in Figure 5, where the reference band value (red 
point) is the fractional rms. We have attempted to obtain 
as many constrained bins as possible, however, the lack of 
correlated variability makes this extremely difficult leaving 
us with only a small number of bins in each observation. 
Irrespective of this we can quite clearly see that the refer- 
ence bandpass has the highest variability (by design), and 
that there is significantly less correlated variability at other 
energies in observations 1, 4 and 5. Additionally, the typi- 
cal energy at which this variability is seen increases as the 
source declines. We note that observation 2 does not neces- 
sarily show a lack of correlated variability as the only other 
well constrained bin of excess variance (besides the reference 
band) is consistent within 1 sigma with the reference band 
whilst the nature of variability in observation 3 is difficult 
to accurately and reliably constrain due to the relative lack 
of counts above 2 keV. 

We accept that it is possible that our iterative energy 
binning scheme could distort the statistical significance of 
this. Indeed, testing this approach on a similarly bright but 
constant source with similar observation length (G21.5-09, 
OBSID: 0122700101) can produce a single narrow bin of 
variability due to the inherent stochasticity of the white 
noise. However, it is highly unlikely that such a process 
could produce the observed energy dependence with lumi- 
nosity soon in M31 ULX-1 and so wc claim that there is real 
variability present. This energy dependence is difficult to 
explain in the single disc model. A disc dominated XMM- 
Newton spectrum from LMC X-3 shows variability on the 
few percent level (Kolehmainen et al. 2011, in prep), but 
this is constant across the entire energy band. Instead, it is 
much easier to explain in the two component model, where 
the interface between the disc and wind picks out a dis- 
tinct radius and hence temperature, and this temperature 
increases as the wind photosphere radius decreases during 
the decline (see Fig 4) . The modulation could be due either 
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Figure 3. Simultaneous X-ray spectral fitting for M31 ULX-1 using the EPIC PN data across all 5 XMM-Newton observations (Obsl: 
black, Obs2: red, Obs3: green, Obs4: blue and Obs5: cyan). Left: we initially attempt to describe the data using a single disc model 
(bhspec). This can provide an acceptable fit quality but only for neutral absorption consistent with the Galactic column to M31. This 
model also suffers from systematically underestimating the flux at soft energies as can be seen in the ratio of data to the model (lower 
panel). This may indicate that the model is too simple, with the residuals potentially indicating emission from within the optically thin 
plunging region. Right: the same data but with a model of slim disc and inner photosphere/ wind production. It is clear from inspecting 
the ratio to the best fitting model that the description is an improvement. This model can also describe the data for a larger absorption 
column (see Table 2). 
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Figure 4. Left: Contribution of the disc component to the best fitting model of slim disk and photosphere/ wind (the total model is given 
as a solid line and the disc component as a dot-dashed line with the colours corresponding to the observations in the caption of Fig. 2). 
As the source dims, the model predicts an increase in advection dominance and a broadening of the component. Right: The contribution 
to the same model from the photosphere/ wind (dot-dashed line with colours as before). As the source dims, the model predicts that this 
component should get hotter and less dominant as it is driven from smaller, hotter radii. 
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Figure 5. Fractional covariance for each of the observations on 300 s binning (to ensure adequate statistics) with the reference band 
shown in red (for which the value is the fractional excess variance). There is clearly a lack of evidence for correlated variability in any 
of the observations, whilst the well-constrained bin of variability is seen to increase in energy. This presents problems for our standard 
view of accretion disc behaviour but can be more readily explained by our ULX model. 



to instabilities in the disc at this radius, or to the wind being 
clumpy and stochastically covering this inner region of the 
disc. 



5 DISCUSSION & CONCLUSION 

The combined spectral and timing behaviour of M31 ULX- 
1 allows a unique insight into the properties of a bona fide 
ULX, albeit a somewhat unusual ULX in being a transient 
rather than persistent system. 

Garcia et al. (2010, ATEL #2474) determined the pres- 
ence of an optical counterpart to the X-ray emission during 
the outburst, on the basis of a ~ 4 ks HST Advanced Cam- 
era for Surveys (ACS) F435W filter exposure taken on 2010 
January 21, with an apparent magnitude of 23.8 (Figure 6, 
left hand panel). This is significantly brighter than a limit 
we derive from a later, post-outburst observation taken on 
2010 July 20 in the same instrument, of mF435w ~ 26 (Fig- 
ure 6, right hand panel). By converting the neutral absorp- 
tion column measured from the X-ray data into an extinc- 
tion using the relation of Predehl & Schmitt (1995), we find 
Ab = 0.66 towards the source in M31 (using the B filter as 
a close proxy of the HST F435W filter). We use this, and 
the magnitude limit, to derive an upper limit on the ab- 
solute magnitude of the optical counterpart in quiescence, 
presumably the emission of the secondary star in the sys- 
tem, of Mf435w ~ +1. This rules out a high mass O or B 
star companion to the compact object, consistent with the 
identification of this transient as an LMXB. The lack of a 
persistent bright optical counterpart also implies that the 
detected transient optical emission is most likely from re- 
processing in the X-ray illuminated outer disc (van Paradijs 
1996). 

The source lightcurve also strongly supports an LMXB 
identification (Fig 2). There is a clear, dramatic outburst, 
followed by an exponential decay. Such outbursts are trig- 
gered by the hydrogen ionisation instability in the outer disc, 
and so requires a low mass transfer rate from the compan- 
ion star. A high mass companion filling its Roche lobe would 
have too large a mass transfer rate, making the system per- 
sistent rather than transient (King 1999). 



The e-fold timescale of 39 days is also very similar to 
that seen in many stellar mass BH LMXB transients in 
our own Galaxy (Chen, Shrader & Livio 1997; Shabaz et 
al. 1998; Mineshige, Yamasaki & Ishizaka 1993). We can 
use this timescale to estimate the radius of the ionised disc 
(from equation 14 of Shahbaz, Charles & King 1998) us- 
ing the peak luminosity estimate from the Chandra observa- 
tion (~5xlO^^ erg s~^), the upper limit for the e-fold time 
(39 days), the accretion efficiency (~0.1 for an alpha disc: 
Shakura & Sunyaev 1973) and the critical density for the 
disc to give an upper limit for the size scale. This gives a 
radius of 2.6x10^^ cm and a lower limit for the period of 
the binary orbit (from Frank, King & Raine 2002) of >32.8 
hours. In terms of binary parameters and inferred luminos- 
ity, this system is then qualitatively similar to 4U 1543-47 
(Chen et al. 1997), reinforcing our identification of the sys- 
tem as a LMXB. For a low mass star to fill its Roche lobe in 
such a wide binary means that the companion is probably 
somewhat nuclear evolved, i.e. a sub-giant. 

The sequence of spectra during the decline can be very 
roughly described by a standard disc model i.e. inner radius 
fixed at the innermost stable circular orbit, around a stellar 
mass (~10 Mq) BH. If this is the correct interpretation then, 
for the first time, we are observing a disc in the Eddington- 
regime without the hindrance of a large neutral Galactic 
column, typical of XRBs in our Galaxy. However, this gives 
residuals of order 10% at low and high energies. This could 
indicate that even our best disc models are incomplete, e.g. 
they do not include advection or release of advected radi- 
ation from the plunging region or self-shielding (Sadowski 
2009; Sadowski et al. 2009). However, these should all give 
an increasing distortion of the disc spectrum at Eddington 
and beyond, yet the residuals are seen at a similar or even 
greater extent in the lowest Eddington fraction spectrum as 
well as in the highest ones. If this is the correct interpre- 
tation of the spectra then our best current disc models can 
only describe the data within 10% over this wide bandpass. 
New disc models that fully incorporate the effects of advec- 
tion, relativistic smearing and energy transport will identify, 
in a more physically robust manner, how current disc mod- 
els fail to describe the emission. This dataset will therefore 



© 2011 RAS, MNRAS 000, 1-10 



The missing link: a LMXB in M31 seen as a ULX 9 



ACS/WFC F435W 
2010/01/21 



1" 



• 



ACS/WFC F435W 
2010/07/20 



1" 



Figure 6. HST (ACS) images taken on the 21st of Jan 2010 (left, OBSID: JB9D15010, PropID=11833) and 20th of July 2010 (right, 
OBSID: JB9D20010, PropID=11833) for exposures of 4360 s each. The highlighted position is that given in the finding charts of Garcia 
et al. (2010, ATEL #2474) and, whilst the former observation corresponds to an apparent optical magnitude of 23.8 we determine the 
latter, quiescent observation to be far fainter: m,p435w ^ 26. Using the Galactic line-of-sight extinction (Predehl & Schmitt 1995), we 
derive an upper limit on the absolute magnitude of the optical counterpart in quiescence of Afp435w ~ +1, ruling out a high mass O or 
B star companion to the compact object. 



prove an invaluable measure of these future models' success 
(Straub et al. in prep). 

Alternatively, the spectra are statistically better fit by 
an advective (slim) disc model, together with a low tem- 
perature, optically thick Comptonisation component which 
could be the photosphere of an inner-disc wind (e.g. Mid- 
dleton et al. 2011b). The observed decrease in extent of this 
wind/photosphere as the luminosity decreases is then due to 
the reduction in mass loss rate into the wind, so that its pho- 
tosphere covers less of the inner disc. The observed energy- 
dependent variability (weakly) supports this two component 
model, but we caution that this is close to the statistical 
limit of what the variability data can probe. 

We note that the situation remains partially degener- 
ate in terms of spectral components and viewing angles. We 
would obtain similar fit statistics if the cool component was 
in fact an outer photosphere and the hot region a modified 
'bare' disc where any loose material has already been up- 
Ufted and removed in a wind. This is our proposed model for 
NGC 5408 X-1 (Middleton et al. 2011a) and could hkewise 
predict the same presence of uncorrelated variability via ex- 
trinsic means but at a different viewing angle. Whilst this is 
a much more likely solution for the higher luminosity ULXs, 
we are currently unable to break this degeneracy in lower 
luminosity ULXs such as this (and M33 X-8; although that 
particular source has no constrained variability: Middleton 
et al. 2011b). In either case, we require a two component 
solution comprising a wind/photosphere and modified disc 
emission. 

Hence we have a transient system that displays an 
extraordinary combination of ULX-like luminosities, an e- 
fold time and derived binary system parameters, and opti- 
cal counterpart consistent with a classical LMXB transient, 
and super-Eddington outburst behaviour similar to other 
ULXs. Additionally, the presence of an IMBH is rejected 
on the basis of the X-ray characteristics, firmly demonstrat- 
ing that an IMBH is not necessary to reach ULX luminosi- 
ties. First and foremost, this then provides a solid eviden- 
tial link between accretion onto an apparent stellar mass 
black hole binary and ULX-like behaviour. This evidence 



exceeds that from previous analyses where the success of 
super-Eddington emission models lead to the inference of 
an underlying stellar-mass black hole (e.g. Gladstone et al. 
2009); in this case the transience allows far stronger physical 
links to be made. 

However, the transient nature of this object also marks 
it out as an unusual ULX. Although it has long been known 
that ULXs are associated with both old and young stellar 
populations (Humphrey et al. 2003; Colbert et al. 2004) , the 
majority of work has focused on the objects directly linked 
to young stellar populations as they are both more numerous 
and more luminous than the older population (cf. Walton et 
al. 2011). These objects are also typically persistently lumi- 
nous, unlike the transient discussed here; in fact few tran- 
sient ULXs have been observed, with notable cases associ- 
ated with the old population of Cen A (Ghosh et al. 2006; 
Sivakoff et al. 2008) and a recurrent transient in the outer 
regions of the archetypal nuclear starburst galaxy NGC 253 
(Bauer & Pietsch 2005). However, despite the likely differ- 
ences in donor star and outburst timescales, the similarity 
in X-ray spectral characteristics between M31 ULX-1 and 
other well-studied low luminosity ULXs (e.g. M33 X-8, Mid- 
dleton et al. 2011b) argues for a commonality in accretion 
physics and an underlying black hole mass in the stellar 
regime in many ULXs. 
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